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The  purpose  of  the  work  presented  in  this  study  is  related  to  heat  transfer  and  airflow  modelling  analysis 
in  solar  chimneys,  according  to  some  dominant  parameters.  A  typical  case  of  application  is  given  in  this 
study.  It  consists  in  analyzing  a  natural  laminar  convective  heat  transfer  problem  taking  place  in  a 
chimney.  Heat  transfer  and  fluid  dynamic  aspects  of  the  airflow,  through  an  axis  symmetric  system  in  a 
dimensionless  form,  with  well  defined  boundary  conditions  is  thus  examined.  Results  are  related  to  the 
temperature  distribution  and  the  velocity  field  in  the  chimney  and  in  the  collector,  determined  by 
solving  the  energy  equation,  and  the  Navier-Stokes  equations,  using  finite  volume  method.  The 
numerical  code  based  on  this  modelling  is  validated  through  the  Vahl  Davis  benchmark  solution  for 
natural  convection  and  to  other  authors  for  other  cases. 

©  2010  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Solar  chimney  power  plant  is  composed  of  a  solar  collector  whose 
function  is  to  increase  the  energy  level  of  the  air  by  greenhouse 
effect,  of  a  chimney  tower  to  ensure  the  air  circulation  through 
density  gradient  and  an  aero  generator  to  produce  electric  power.  It 
is  a  natural  generator  of  power  which  uses  solar  radiation  to  increase 
the  internal  energy  of  the  air  flowing  in  the  system,  and  then  to 
transform  the  useful  profit  of  the  solar  collector  into  kinetic  energy  of 
flow  which  might  be  transformed  into  electric  power  by  means  of  a 
suitable  turbine.  Since  the  first  investigations  of  Schlaich  [1]  and 
Schlaich  et  al.  [2]  in  the  solar  chimney  field,  efforts  were  focused 
essentially  on  analyzing  performances  and  cost  of  solar  chimney 
power  plants  (SCPP).  It  should  be  noted  that  the  available  literature 
on  fluid  dynamics  aspect  analysis,  and  on  the  one  related  to 
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transport  equations  solution,  given  by  computational  fluid  dynam¬ 
ics  is  scares  because  investigations  have  concentrated,  most  of  the 
time,  on  the  total  energy  performances  and  cost  evaluation  of  such 
systems  [3-14].  Fluid  dynamic  aspects  of  these  problems  related  to 
the  local  characteristic  study  of  the  flow  become  better  known,  such 
as  detection  of  spades  and  recirculations  zones  as  well  as  weak 
temperature  variations.  Bernardes  et  al.  [  1 5  ]  conducted  a  theoretical 
analysis  of  a  solar  radial  air  heater,  operating  on  natural  laminar 
convection  in  steady  state,  to  predict  the  thermo-hydrodynamic 
behaviour  of  the  solar  chimney  device.  Finite  volume  method  in 
generalized  coordinates  was  used  to  solve  the  governing  equations. 
Velocity  field  and  temperature  distribution  in  the  flow  were 
obtained  under  imposed  thermal  conditions.  Von  Backstrom  and 
Gannon  [16]  were  interested  mainly  in  a  one-dimensional 
compressible  flow  for  the  thermodynamic  variable  calculation 
versus  the  tower  height.  Von  Backstrom  and  Fluri  [17]  investigated 
analytically  the  validity  and  applicability  of  the  assumption  that,  for 
maximum  fluid  power,  the  optimum  ratio  of  turbine  pressure  drop 
to  pressure  potential  is  2/3.  Pastohr  et  al.  [18]  used  the  FLUENT 
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Nomenclature 

cp  specific  heat  (J/kg  K) 

D  diameter  (m) 

e  height  of  the  cover  (m) 

Gr  Grashoff  number 

h  heat  transfer  coefficient  (W/m2  K) 

H  height  (m) 

Nu  Nusselt  number 

p  pressure  (Pa) 

P  power  (W) 

Pr  Prandt  number 

Ra  Rayleigh  number 

r,  y  space  coordinates  (m) 

t  time  (s) 

T  temperature  (K) 

u,v  velocity  components  (m/s) 

Greek  symbols 

X  thermal  conductivity  (W/m  K) 

p  density  of  air  (kg/m3) 

pt  viscosity  of  air  (kg/m  s) 

/3  thermal  coefficient  of  expansion  (K_1) 

Indexes  and  subscripts 
h  hot 

c  cold 

cov  cover 

soil  soil 

0  inlet 

ref  reference 


software  for  modelling  a  solar  chimney  power  plant,  geometrically 
similar  to  that  of  Mansaranes,  with  the  aim  of  carrying  out  an 
analysis  and  reporting  details  on  the  operating  mode  and  the  system 
efficiency.  They  confirmed  that  the  pressure  drop  in  the  turbine  and 
the  mass  flow  rate,  decisive  elements  on  the  system  effectiveness, 
cannot  be  only  given  by  coupling  all  the  parts  of  a  SCPP.  Numerical 
results  given  by  FLUENT  software  are  in  good  agreement  with  the 
results  given  by  a  simple  model  proposed  by  the  authors,  which  led 
to  the  conclusion  that  it  is  much  easier  to  use  it  for  parametric 
studies.  Pretorius  and  Kroger  [19]  evaluated  the  influence  of  a 
developed  convective  heat  transfer  equation,  more  accurate  turbine 
inlet  loss  coefficient,  quality  collector  roof  glass  and  various  types  of 
soil  on  the  performance  of  a  large  scale  solar  chimney  power  plant. 
Ming  et  al.  [20]  presented  a  mathematical  model  to  evaluate  the 
relative  static  pressure  and  driving  force  of  the  solar  chimney  power 
plant  system  and  verified  the  model  with  numerical  simulations. 
Later,  Ming  et  al.  [21  ]  developed  a  comprehensive  model  to  evaluate 
the  performance  of  a  solar  chimney  power  plant  system  in  which  the 
effects  of  various  parameters  on  the  relative  static  pressure,  driving 
force,  power  output  and  efficiency  have  been  further  investigated. 
The  authors  supposed  the  existing  models  are  insufficient  to 
accurately  describe  all  the  phenomena  occurring  in  solar  chimney 
power  plant.  Using  the  solar  chimney  prototype  of  Manzanares,  as  a 
practical  example,  3-D  turbulent  flow  numerical  simulation  studies 
were  performed  to  explore  the  geometrical  modifications  on  the 
system  performance.  Results  showed  a  good  agreement  with  the 
analytical  model.  The  control  of  the  SCPP  analytical  tools  such  as 
dynamic  simulation  of  these  systems  is  essential.  Numerical 
simulations  were  performed  by  Ming  et  al.  [22]  to  analyze  the 


characteristics  of  heat  transfer  and  airflow  in  the  solar  chimney 
power  plant  system  with  energy  storage  layer.  Different  mathemat¬ 
ical  models  for  the  collector,  the  chimney  and  the  energy  storage 
layer  were  established,  and  the  effect  of  solar  radiation  on  the  heat 
storage  characteristic  of  the  energy  storage  layer  was  analyzed. 
Numerical  simulation  results  show  that  the  heat  storage  decreases 
firstly  and  then  increases  with  the  increase  of  the  solar  radiation 
from  200  W/m2  to  800  W/m2.  The  static  pressure  decreases  while 
the  velocity  increases  significantly  inside  the  system  with  the 
increase  of  solar  radiation;  the  average  temperature  at  the  outlet  of 
the  chimney  and  the  one  of  the  energy  storage  layer  may  increase  too 
significantly  with  the  increase  in  solar  radiation.  In  addition,  the 
temperature  gradient  of  the  storage  medium  may  increase  and  this 
results  in  an  increase  of  energy  loss  from  the  bottom  of  the  energy 
storage  layer.  Pastohr  et  al.  [23]  presented  a  numerical  simulation 
result  in  which  the  storage  layer  was  regarded  as  solid.  In  their  paper, 
conjugate  numerical  simulations  of  the  energy  storage  layer,  the 
collector  and  the  chimney  have  been  conducted,  and  the  character¬ 
istics  of  the  heat  storage  system,  and  the  flow  and  heat  transfer  in  the 
whole  system  have  been  studied.  Zhou  et  al.  [24]  have  performed  an 
experimental  study  in  a  solar  chimney.  A  pilot  experimental  solar 
chimney  power  setup  consisting  of  an  air  collector  of  10  m  in 
diameter  and  an  8  m  tall  chimney  was  built.  The  authors  noted  that 
the  temperature  difference  between  the  collector  outlet  tempera¬ 
ture  and  the  one  of  the  ambient  usually  might  reach  as  much  as 
24.1  °C,  which  generates  the  driving  force  of  the  airflow  in  the  setup. 
Their  data  analysis  showed  an  air  temperature  inversion  in  the  latter 
chimney  after  sunrise  and  this  is  due  to  the  increase  of  solar 
radiation  from  the  minimum.  The  phenomenon  clears  up  once  a 
driving  force  is  generated  by  a  temperature  high  enough  to 
overcome  it.  Maia  et  al.  [25]  presented  a  theoretical  analysis  of  a 
turbulent  flow  inside  a  solar  chimney.  They  showed  that  the  most 
important  physical  elements  in  a  solar  chimney  system  are  the 
tower  dimensions  as  they  cause  the  most  significant  variation  in  the 
flow  behaviour.  An  increase  in  the  height  and  in  the  diameter  of  the 
tower  produces  an  increase  in  the  mass  flow  rate  and  a  decrease  in 
the  flow  temperature.  Chergui  et  al.  [26,27]  simulated  a  thermo¬ 
hydrodynamic  behaviour  analysis  of  the  airflow  through  an 
axisymmetric  system,  such  as  chimneys,  with  defined  boundary 
conditions.  Special  emphasis  is  given  to  an  example  of  application  of 
natural  laminar  convective  heat  transfer  problem  occurring  in  a 
solar  chimney  power  plant. 

The  above  work  focused  mainly  on  global  analyses  which  were 
developed  on  this  type  of  system  include  calculations  of  output 
energy,  system  efficiency,  parametric  analyses  and  analytical 
models.  Mostly  the  results  showed  a  viable  system  technically  and 
economically  in  using  solar  energy  resources.  However  this 
literature  review  showed  a  lack  of  information  on  the  local 
characteristics  of  the  flow  inside  the  collector-chimney  system, 
information  needed  for  system  configuration  and  design  in  order  to 
in  consideration  mechanical  conversion  system. 

In  order  to  bring  some  more  understanding  to  the  phenomenon 
taking  place  during  the  flow  through  the  collector  and  the  chimney. 

Using  the  finite  volume  method  under  certain  assumptions,  the 
velocity  field  and  temperature  distribution  are  determined  by 
solving  the  mathematical  model  represented  by  the  Navier-Stokes 
equations,  continuity  and  energy  equations,  in  dimensionless  form 
which  will  give  us  local  characteristics  of  the  flow  and  let  us 
determine  the  suitable  zone  for  the  positioning  of  mechanical 
conversion  system. 

2.  Governing  equations  and  boundary  conditions 

The  configuration  model  represents  a  cylindrical  cavity  with  a 
ring-like  inlet  opening  at  the  bottom  and  a  circular  outlet  opening 
at  the  top.  Assumptions  are  constant  properties  except  in  the 
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formulation  for  the  buoyancy  term  in  which  Boussinesq  approx¬ 
imations  were  made,  negligible  compressibility  effects  and  viscous 
dissipation,  Newtonian  fluid,  and  laminar  and  two-dimensional 
flow.  We  assume  an  axisymmetric  flow. 

The  governing  equations  are  given  by: 


-  Continuity  equation: 

|(P)+}|(P™)+|(pv)  =  o 


Momentum  equations: 

9  /  X 
(Pu)  + 

dp 


d  19  9 

di(pu)+7dr  (pr“2)  +  dy (puv) 


id/  du 
dr  ^  r  drx^  dr 


9 


du 


+  dy  V  3 Y. 


U 

2p-j 


(1) 


(2) 


i(pv)  +  7ppruv)+^(pv2) 

dp  19/  dv\  9  (  dv\  .  . 

=  ”9y+7  dr{rfldr)  +3y(/3j/  +{p0~p)g 


-  Energy  equation: 

9  1  9  9  1  d  f  X  dT 

9t  (pT)  +  7  3F (pruT)  +  9y (pl/T)  =  7  ( r  Z  97 


(3) 


.  2.  (2  HE 

+  dy  \cp  dy. 


(4) 


The  use  of  the  dimensionless  variables  for  the  previous 
equations  allows  the  generalization  of  the  result  which  can  be 
extended  to  any  dimensional  system. 

Let  us  then  consider  H  a  linear  characteristic  dimension  of  the 
flow,  uref  a  reference  velocity,  and  AT  a  reference  temperature 
difference. 

The  dimensionless  variables  will  then  be  given  by: 


t+  = 


t  +  x  +  y  + 

—  ,x+=-,y+=fj,u+ 

Lj*0^  il  fi 


U  ,  V 

,V+  = 


’+ 


llref 


Ziref 


T-T 


A  T 


ref ,  p+ 


P  ref 


In  the  free  convective  heat  transfer  case,  the  reference 
parameters  are  defined  by  [4]: 

,  H  /  Oj_J  \  T  rji  (Tq  +  Th) 

fref  =  5  =  V^HAT’  =  2  ’ 

=  pulf  (6) 

with  AT  =  Th  -  Tc. 

Thus,  the  dimensionless  variables  become  then: 


i-+ 


t\i^2,x+  =^,y+ 


H’  /iGr 


UPH  ,,+ 

0.5  ’  V 


VpH 

/xGr1 


0.5 


(7) 


T+  =  7  ,/ef,p+  =  1 Gr  =  pr  =  fp- , Pa 


AT  ’  p.2Gr 
p2g/3ATH3Pr 

P'2 


PL' 


X 


(8) 


The  non-dimensional  general  form  of  the  differential  equation 
describing  the  flow  is  given  by: 


1  3  7  .  r0+  #_ 

r+  9r+  \  9r+ 


+  x  9 
+ 


9y 


+ 


r 


0+  90 


+ 


dy 


+ 


+  S1 


(9) 


Table  1 

Values  of  @+,  r(p+  and  S0+  of  Eq.  (9). 

Equation  @+ 

Mass  balance  / 

Momentum  balance  versus  r  component  u+ 

Momentum  balance  versus  y  component  v+ 

Energy  T* 


T57 

o 


dp+ 

dr+ 

d  p+ 
dy+ 


2u+ 

~  r^~ 

+  T+ 


0 


5u/  5y  =  0 
5v/  5y  =  0 
5T/6y=0 


0  stands  for  1,  u,  v  or  T  according  to  the  conservative  governing 
equation  considered  and  r*  and  S ^  are  the  corresponding  diffusion 
coefficient  and  source  term  respectively  (Table  1 ). 

However,  the  study  domain  is  delimited,  in  the  cylinder,  in  using 
a  numerical  subtlety  consisting  in  taking  a  very  high  flow  viscosity  in 
the  volume  outside  the  collector-chimney  system  (Fig.  1). 

The  solution  to  the  generated  set  of  equations  must  satisfy  some 
conditions.  Fig.  1  shows  the  physical  model.  Tcov  is  the  cover  tem¬ 
perature.  The  chimney  wall  is  considered  adiabatic  ((9T/9n)|w  =  0). 
The  non-slip  condition  is  imposed  on  all  walls  (u  =  v  =  0).  TSOii  is  the 
ground  surface  temperature.  The  cover  and  the  ground  tempera¬ 
tures  are  obtained  by  energy  balance  [25].  Since  an  axisymmetric 
flow  is  assumed  in  the  chimney,  then  u  =  0,  (dv/dr)  =  0,  (dT/dr)  =  0 
at  r  =  0.  At  the  top  of  the  chimney  tower,  a  fully  developed  flow 
region  is  assumed  on  the  velocity  and  on  the  temperature.  At  the 
inflow  boundary,  the  entrance  temperature  T0  is  assumed  constant. 
As  it  is  specified  by  Bernardes  et  al.  [15],  in  natural  convection 
problems,  the  mass  flow  at  the  entrance  caused  by  buoyancy  forces 
is  unknown  in  the  beginning.  The  u-component  of  the  velocity  at  the 
entrance  is  updated  at  each  iteration  by  the  values  of  the 
neighboring  velocities.  For  the  initial  conditions,  u  =  0;  v  =  0;  T,  =  0. 

3.  Numerical  method 

A  finite  volume  method  [28,29]  is  used  to  discretize  the 
governing  equations.  The  calculation  domain  is  uniformly  divided 
into  small  non-overlapping  control  volumes.  A  staggered  grid  is  used 
such  as  the  velocities  lie  on  the  faces  of  the  control  volumes  whereas 
the  pressures  and  the  temperatures  are  located  at  the  centers.  A 
primitive  variable  formulation  is  used  to  solve  the  equations  and  a 
fully  implicit  time  is  employed.  The  power  law  scheme  is  used  to 
discretize  the  convection  and  diffusion  terms  and  the  SIMPLE 
algorithm  is  used  to  solve  the  problem. 
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Fig.  4.  Dimensionless  isovelocity  lines  for:  Ra  =  106  and  e/H  =  0.1. 


The  worked  out  computer  program  is  validated  with  respect  to 
the  de  Vahl  Davis  benchmark  solution  [30]  and  others  [31,32].  Of 
course,  the  generated  program  was  validated  over  a  flow  in  a 
cylindrical  cavity  for  different  flow  regime  [33]. 

The  grid  dependence  has  been  investigated  using  different 
mesh  sizes  before  settling  to  a  mesh  size  of  1/22  (24  x  24  cells), 
although  the  computer  program  was  validated  using  mesh  sizes 
identical  to  those  used  in  the  reference  work. 

As  far  as  the  time  step  is  concerned,  different  values  of  At  were 
considered.  Slight  variations  of  the  calculated  functions  appear 


around  t=  1.  However  the  effect  of  the  time  step  size  becomes 
negligible  when  At  <  0.2  (the  deviation  of  the  average  temperature 
appears  on  a  small  zone  around  t  =  1 .  and  is,  at  most,  equal  to  1 .4%). 

Since  we  consider  the  influence  of  the  mass  flow  rate  through 
the  variation  of  the  Reynolds  number,  a  test  was  made  as  far  as  the 
flowing  regime  is  concerned.  In  this  case,  we  took  Re  =  2.3  x  10+4. 

The  simulation  is  made  until  the  steady  state  is  reached  then  we 
analyzed  the  evolution  of  the  velocity  components.  It  appears  that, 
once  the  steady  state  is  gotten,  both  components  became  constant 
confirming  the  assumption  of  laminar  flow.  The  validation  was 
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Fig.  3.  Dimensionless  isovelocity  lines  for:  Ra  =  104  and  e/H  =  0.1. 


Fig.  5.  Dimensionless  isovelocity  lines  for:  Ra  =  108  and  e/H  =  0.1. 


1414 


T.  Chergui  et  al. / Renewable  and  Sustainable  Energy  Reviews  14  (2010)  1410-1418 


Fig.  6.  Dimensionless  isothermal  lines  for:  Ra  =  100  and  e/H  =  0.1. 


Fig.  8.  Dimensionless  isothermal  lines  for:  Ra  =  106  and  e/H  =  0.1. 


made  using  water  as  working  fluid  flowing  into  a  cylindrical  cavity. 
For  more  details  concerning  this  question,  the  reader  is  referred  to 
the  study  in  [33]. 

4.  Results  and  discussion 

Since  the  governing  equations  are  written  in  non-dimensional 
form,  the  influence  of  the  fluid  properties  and  the  flow  rate  are 
studied  through  the  Prandtl  number  and  the  variation  of  the 
Rayleigh  number.  We  consider  e/H  =  1.  The  case  study  consists  in 


analyzing  a  natural  laminar  convective  heat  transfer  problem 
taking  place  in  a  chimney.  Heat  transfer  and  fluid  dynamic 
aspects  of  the  airflow,  through  an  axis  symmetric  system  in  a 
dimensionless  form,  with  boundary  conditions  defined  in  Section 
2  is  thus  examined.  The  results  are  related  to  the  temperature 
distribution  and  the  velocity  field  in  the  chimney  and  the 
collector. 

Figs.  2-5  illustrate  the  dimensionless  velocity  lines  for  e/ 
H  =  0.1  and  for  Rayleigh  numbers  equal  to  100,  104,  106  and  108 
respectively.  It  should  be  noted  that  the  velocity  magnitude 


Fig.  7.  Dimensionless  isothermal  lines  for:  Ra  =  104  and  e/H  =  0.1. 


Fig.  9.  Dimensionless  isothermal  lines  for:  Ra  =  108  and  e/H  =  0.1. 
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Fig.  10.  Dimensionless  isovelocity  lines  for:  Ra  =  100  and  e/H  =  0.01. 


Fig.  13.  Dimensionless  isovelocity  lines  for:  Ra  =  108  and  e/H  =  0.01. 
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Fig.  11.  Dimensionless  isovelocity  lines  for:  Ra  =  104  and  e/H  =  0.01. 


Fig.  14.  Dimensionless  isothermal  lines  for:  Ra  =  100  and  e/H  =  0.01. 
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Fig.  12.  Dimensionless  isovelocity  lines  for:  Ra  =  106  and  e/H  =  0.01. 


Fig.  15.  Dimensionless  isothermal  lines  for:  Ra  =  104  and  e/H  =  0.01. 
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Fig.  16.  Dimensionless  isothermal  lines  for:  Ra  =  106  and  e/H  =  0.01. 


Fig.  18.  Velocity  distribution  versus  collector  radius  for:  Ra  =  100  and  e/H  =  0.01. 


increases  with  the  increase  of  the  Rayleigh  number  and  its 
maximum  value  is  located  approximately  at  the  inlet  of  the 
chimney  as  it  is  reported  in  some  literatures  [  1 8,21  ].  For  all  these 
Rayleigh  numbers,  velocity  lines  are  smooth  characterizing  a 
laminar  flow,  along  the  chimney  and  the  collector,  except  for  the 
Rayleigh  number  of  108  where  different  zones  of  flow  recircula¬ 
tion  can  be  observed  in  the  collector  as  well  as  in  the  chimney  due 
to  instability  of  airflow  in  these  areas.  These  recirculation  zones 
remind  us  the  Rayleigh-Benard  flows  beyond  a  critical  value  of 
Rayleigh  numbers. 

Figs.  6-9  show  the  dimensionless  isothermal  lines  for  e/H  =  0.1 
and  for  Rayleigh  numbers  equal  to  100,  104,  106  and  108 
respectively.  The  maximum  temperature  is  located  near  the 
ground  in  the  collector  due  to  heat  transfer  exchange  between  this 
surface  and  the  airflow  beneath  the  cover.  For  a  Rayleigh  number 
of  108  swirls  in  isothermal  lines  is  also  observed  in  the  collector 
area. 
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Figs.  10-13  show  the  dimensionless  velocity  lines  for  e/H  =  0.01 
and  for  Rayleigh  numbers  equal  to  100,  104,  106  and  108 
respectively.  Notice  that  the  velocity  magnitude  increases  with 
Rayleigh  numbers  and  its  maximum  is  located  approximately  at 
the  inlet  of  the  chimney.  However  its  value  is  smaller  than  the  one 
calculated  in  the  previous  case  with  e/H  =  0.1.  Velocity  lines  are 
smooth  along  the  collector  and  the  chimney  for  all  Rayleigh 
numbers. 

Figs.  14-17  show  the  dimensionless  isothermal  lines  for  e\ 
H  =  0.01  and  for  Rayleigh  numbers  equal  to  100,  104,  106  and  108 
respectively.  The  maximum  temperature  is  still  located  at  the 
lower  part  of  the  collector.  For  the  Rayleigh  number  of  108,  a  slight 
disturbance  can  be  observed  in  the  collector  zone. 

Figs.  18-21  give  the  dimensionless  air  velocity  versus 
the  dimensionless  collector  radius  for  e/H  =  0.01  and  for 
Rayleigh  numbers  equal  to  100,  104,  106  and  108  respectively. 
We  notice  that,  on  all  these  figures,  the  air  velocity  decreases 
with  the  collector  radius  increase  and  this  magnitude  increases 
with  Rayleigh  number.  The  same  results  are  observed  by  Pastohr 


Fig.  17.  Dimensionless  isothermal  lines  for:  Ra  =  108  and  e/H  =  0.01. 


Fig.  19.  Velocity  distribution  versus  collector  radius  for:  Ra  =  104  and  e/H  =  0.01. 
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Fig.  20.  Velocity  distribution  versus  collector  radius  for:  Ra  =  106  and  e/H  =  0.01. 


Dimensionless  Collector  Radius 


Fig.  21.  Velocity  distribution  versus  collector  radius  for:  Ra  =  108  and  e/H  =  0.01. 


et  al.  [18]  in  their  study  in  which  dimensional  equations  are 
used. 

5.  Conclusion 

Many  encouraging  studies  were  reported  about  the  use  of  solar 
chimney  to  produce  electrical  power.  Performance  studies  were 
related  to  geometrical  and  operating  parameters.  This  study 
considers  the  heat  transfer  process  and  the  fluid  flow  in  the 
collector  and  the  chimney  under  some  imposed  operational 
conditions.  A  validated  computer  programme  was  adapted  to  the 
solar  chimney  configuration  to  solve  the  primitive  governing 
equations  (continuity,  momentum  and  energy  equations).  Con¬ 
sidering  the  collector  height  and  the  driving  force  of  the  flow  i.e. 
the  temperature  difference  through  the  Rayleigh  number,  the 
velocity  field  and  the  temperature  distribution,  through  the 
system,  are  evaluated.  Results  showed  local  flow  characteristics 
and  it  appears  that  for  most  Rayleigh  numbers,  the  flow  seems 
laminar  except  for  Raleigh  number  of  108  where  there  are  some 
disturbances. 


Maximum  velocities  are  gotten  at  the  entrance  of  the  chimney 
tower  and  its  magnitude  is  higher  in  the  case  of  e/H  =  0.1  for  all 
Rayleigh  numbers. 

A  follow  up  of  the  velocity  field  in  the  collector  shows  an 
increasing  magnitude  from  the  periphery  to  the  center.  This  result 
will  let  the  solar  chimney  designer  correctly  locate  the  mechanical 
conversion  device  inside  the  collector  chimney  system. 
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